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THE OSC ILLA TING CIRCULAR AIRFOIL ON THE BASIS OF POTENTIAL THEORY 1 

PART I t 
By Th. Schade 


Proceeding from the thesis hy W. Kinner the present report treats 
the prohlem of the circular airfoil in uniform airflow executing small 
oscillations, the amplitudes of which correspond to whole functions of 
the second degree in x and y. The pressure distribution is secured 
hy means of Prandtl's acceleration potential. It results in a system 
of linear equations the coefficients of which can he calculated exactly 
with the aid of exponential functions and Hankel's functions. The 
equations necessary are derived in part I; the numerical calculation 
follows in part II. 


INTRODUCTION 


The present study was undertaken in connection with Kinner' 8 thesis 
(reference l). It deals with the corresponding prohlem of the oscillating 
circular airfoil. 

The specific quantity necessary for predicting the forces and moments 
of an oscillating and deforming airfoil is the z coordinate of the lift- 
ing surface, z = z(x, y, t) . The downwash function W(x, y, t) stands 
with it in the relationship 


dz 

— + 

dt 



W 


Since all oscillations and deformations of the second degree in x and 
y are computed, the total downwash functions necessary thereto are 
represented hy 


1 "Theorie der schwingenden kreisfcSrmigen Tragfl'ache auf potential-- 
theoretischer Grundlage." Luftfahrtforschung, vol. 17, no. 11-12, Dec. 
10, 1940, pp, 387-400, and vol. 19, no. 8, Aug. 20, 1942, pp. 282-291. 
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W = [A + Bx + Cy + Dx 2 + E xy + Fy 2 ] e iVt 

whereby the constants A . . . F may assume any complex values. 

For attire viat ion of the subsequent formulas the following symbolic 
style of writing is introduced 
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+ coy 
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! dz 
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THE VELOCITY POTENTIAL (reference 2) 

The flow velocity of the circular . airfoil is constant .. (V), the 
x-axis in direction of V points rearward, the y— axis in spar direc- 
tion to the right, the z-axis at right angles upward (fig. l) . 

With w denoting the velocity vector, b the acceleration vector 
and t the time, the relation 


b 





is valid. 

Then the introduction of b = grad cp and w = grad <p followed by 
integration yields 


q> = £$' + V 
cJt dx 


As shown later on, cp indicates, up to a certain factor, the pres- 
sing at the upper and lower surface of the disk direct hence cp may also 
be termed the pressure potential. 


Then if the pressure potential cp(x, t) is given, the velocity 


potential $ can be computed from the above equation: 


$ 


1 


(x, y, z, t) ■ i 





t _ 


X - X 


’) 


dx* 


The integration is carried out after putting 


x* = x - V (t - t*) 


where x, y, 0, and t indicate the coordinates of • the starting point. 
The pressure potential cp satisfies Laplace*s equation 


Acp = 0 


whence, after putting 


cp = C (x, y, z) f(t) 
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\|/ (x, y, z) must satisfy the condition A\J/ = 0. For the present f(t) 
is any time function and C a constant. The time function follows from 
the required form of downwash. But since this shall he harmonic in t, 
f (t) itself must he a harmonic function. Then the vertical component of 
the velocity is 



The form of the lifting surface follows from 


v ~ 


+ Y — 
3t cte: 


at 


Z ss 


z 

I W 


■)o 


y, 0, t - 


~) dx» + J?(, y, t - -) 

Y J \ lJ 


where function F serves to assume any position of the wing leading 
edge . 


The pressure potential cp to he defined then must he a discontinu- 
ous function on the circular surface, hut a continuous function at all 
other points in space. 


THE POTENTIAL FUNCTIONS OF THE FIRST KIND 


For the sake of brevity the reader is referred to Earner's report 
(reference 1) where these problems are treated in detail. 

Laplace's equation 


= K 0 

3x 2 dj 2 5z 2 

is transformed to elliptic coordinates through the orthogonal trans~ 
formation 
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x = */ 1 — \/l + q 2 c,os 

y = <7 1 - n a */l + Tp Bin $ 
z = k q 


The x, y, z coordinates are made dimensionless with the "base circle 
radius c. 


With the application of a separation theorem, the real potential 
functions of the first kind with discontinuity at the disk then are 
according to Kinner (reference l): 


m 


* ) « 


.n - m + 1 


<a£g>i eg ig m c(i,)f oe “-n 

. . . 11 n 11 Isinmdj 

(n + m) I •' 


n + m "being posted odd, C™ denoting a constant factor, and factor 

* merely serving to simplify the calculation. 

(n + m) I 


DOWNWASH 


The downwash is then computed for the potential functions previously 
enumerated, resulting in the following downwash expressions: 


W - f (x, T ) e lvt + e l(ut - “> g(,l (3.1) 

where f(x, y) is a whole rational function, hut e ia)X g(y) is trans- 
cendental in the second term. 

The downwash functions are so selected and linearly superimposed 
that f(x, y) exactly represents the l-equired downwash. 


The term containing the factor e u x is then inadequate still 
for representing the disk motion. Subsequently, it is linearly superposed 
with the downwash functions of the second type, so that this total term 
disappears. The solution of this problem is hereinafter termed "com- 
pensation calculus," a detailed discussion of which is given in the sec- 
tion, Compensation Calculus. It affords the linear equation systems, 
the solution of which gives the coefficients for the prediction of the 
pressure distribution. 
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To compute the down-wash, it is necessary to form 


x 



z -> 


— 00 


0 


dx 


(3.2) 


the integration taking place on a streamline. 


With. x Q = ± •J 1 - y 2 denoting the x coordinate of the disk edge and 

^ _ y x 2 + -y 2 + z 2 the radius of any sphere about the disk center, the 
c 

total integration space can be divided in the outer space with r > c 
and the inner space with r < c . 


It proved advantageous to carry out the integration with respect to 
x according to (3.2) separately for the outer and inner space. 

Equation (3*2) is written: 


m 


W 

n 




CO) 

K z 'i — 3> 0 


dx 


m 


= W 


m 


n, a 


+ W 


n.i 


— co 


-Xq 


For the integration over x in the plane of the disk (z = 0) the 
following formulas should be considered: 


1. For the outer space: 


x > 0 


dx 

dp 

dz 


rj d T) 

"/x§ + ti 2 

1 

n 


p = 0 




y 


2 


cos-d = 



" . 0 
dz 
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2. For the inner space: 


dx = T 


/ 




(id4 


2 2 

X c ~4 


n = 0 


cos.£ = 


/*§ - 4 s 


1 - 4 


84 

Sz 


= 0 


Srj 1 
Sz 4 


Sri 

Sz 


= 0 


HEPBESENTATION OF FLAPPING MOTION 


The requisite downvash formula for the flapping motion reads 


W » A e : 


ii’t 


(4.1) 


the pressure potential 


•Pi = C° 4(1 - h arc ctgTi)e 


irn; 


"being used for the representation. 
Then the downwash follows at 


x 


~ o 

0 Cl n Scpx 
Wi = — / lim — dx' 

V / z > o 


which gives for Inner space 


x 


w!,i - 2i P-- e ivt ' dx' 


■X r 


(4.2) 
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00 

and withjv ( t—i ') = co (z — x ')\ from a> = -y- follows 1 Ff.. = J* (l — *. V**— 1 arc ctg z, V**— 1) yp== d» 

x 1 




Wj.— 


_» CJ 1 trt , _«»(»(—“*) 

* . .• & T o ' xr * 


~ CO ixt , 


2 F co i c ^2 F co i 

The dovnvash for the outer space is: 

Wi .. = ^e‘ <vt-<ux) f(l — arc ctg jj) — yy==f^- <*»i- 

For the solution of the integral 
the following substitution is made, and 
this is equally applied on all further 
outer space integrals: 

t)-— x 0 Vs* — 1 V *0 +*?* = x 0 s 
d rj d » 

ui+o' ~ y^r ’ 

In the subsequent calculations the 
C£ i(»fc_cux) 


and, after partial integration: 

if . co 

Wi. • = [*• • are ctg x, f. i*^T- e ~ " * *• *] + 


(4,3) 


ds + —. P 

f s ._l n 01 J + i — 

whence the down wash of the outer space 

l°>-Xo)+ ■ 

1 1 f x < ,s.e- u ’ ix ‘ t d s 

<oi ^ zj ** + 1 xl V 8* — 1‘ 

_ ? V— 

where - y i Hf («, . z.) = J ± =^_ . ds 

denotes Hanked* s cylindrical function. 

The solution of the determinate integral 
over s is achieved by transformation in 
an indeterminate integral over to. Putting 
y = sin a, Xq = cos a, then yields 


*•* 


— aix.t 


dt. 


factor 


is then omitted. 


and indicated by dash over W. 
Herewith: 


? xts.e-”** 1 dt / l 1 \ 

^ xgs 2 -fl— zjj y g > — 1 2jJ 1 / s 1 — 1 \s • cos a -j- i sin a s-cosa — isinaj' 

1 T C e — “* (« cosa-)-i »lna) S? — iui(«-cosa — isina) 

— — e — “ Bln a l — ; d ,-f- f . a,si0a ( - rf - | 

2 L jj (*• cos a + isina) )«* — 1 J (s-cosa — i sin a) y s' — X J 


In this case, however, the relations 


d 

dai 


J e ~ m i (s cos c 
(s • cos a + i si 


03 i (scosa + t sin a) 


sin a) y s 2 — 1 


ds — 


00 

f 


~coi(t COS a -f - 1 sin a) 

y s »^r 


ds = -?- e “ aiDa ff!*>( a ,x t ), 


i sin a) 


j /» - — f*> i (s COS a — isina) r* — 10 i (b • COS a — is , _ 

— T— . . - w— = -< = -f *~“ 8ina (<ox °> - 

da> j (t • cos a — 1 sin a) p s 2 — 1 J p »• — 1 « 

hold true, and with it the above integral becomes: 

f ° x 0 • s • e~ " it|< ds 

3 zfs* + l-zj 'y'iYTT 
1 

CO ox 

= i e _ “ 8in “ [c* 11 — J*| ^ z) («J cos «) - e+ fT!sta a d ToJ + -i- e+ " stn a J y ^ 0 8) (co cosa) •«- s « lna ico] = 

0 0 
to 

= y ^i 1 ’ - e — " sln 0 + y C' 11 - e+ “ ®I“ “ — y J (co • cos a) • CTof [(a, — ai) sin a] d co = 

0 

QJ 

== i Qq (y) ©in (co sin <x) -j- P 0 (y) • (co sin a) — ^o 2) (® cos a ) * ^°f [(^ — a>) sin oc] c? co , (4, 6) 

o 

the constants of integration and c! 1 ^ being computed according 

to (reference 3) • 


d 2K+1, = 

2 


■cf K) = 
2 


J ( S cosa±/sm«)^ + i y 5 T=T _( “ 1)K (2 P « T iQ * K ) 

j (s'cosaiW*)** ys*-l = ( “ 1 ) Jr (<?«- 1±, 7 J, tJr- 1 )J 


(4,6) 


x tD here denotes the real reduced frequency, while a) in the subse- 
quent calculations are integration variables. 
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The total dovnwash la 

and herewith 

"t — i-T •^+9‘* (W -“ X, [-I‘ fl< o«( a . co s«)+ = l_i_ F o | 1 ds 

. L W *-‘ 2id 0 ,^--+2 i JxS^+ 1-xB V7*-T 

-f- J* A?**’ (5> cos a) • <£o[ [(e> — «5) sin a] d tu + 


By comparing this result with the 
outer space integral of the flapping 
oscillation, it can also he written in 
the form 


vhlch means: 


. 1 , . , , 1 4- sin a 

4-j — -®tn (to sin a) -In -r— 1 — : 

'a® ' ' 1 — sin a 


»n=- 


cl 3 


2a> i 


- • ttof (<u ■ sin 


in*)] 


(4,7) 


I ^2 J> (vt — (o 
+ -y~* 


T‘‘-‘ (fj+^) + 


+ 


+ { 


/I 3 

3 d 

f ds 

\,2 '2 to* 

2(0 d to 

,) ** + 1 — *o y*‘— x J 


The indeterminate integral over to 
is definitely computable and is given 
later as series in cd. The first term gives the requisite 

REPRESENTATION OF TORSIONAL OSCILLATION dovnwash and. In addition, a flapping 
ABOUT THE Y-AXIS oscillation, which can he compensated 

. . again hy linear superposition with Wi°, 

(RELAXATION OSCILLATION) Putti ng for this purpose 

The requisite dovnwash formula shall read: „ „ •> t-i 

and forming if i + Tf° = w^. 


W = B■x■e l1 ' , (6,1) 2 

This furnishes the potential function gives the dovnwash 

• • — V X UT W 1 -t— TJ“ A W 

V 4 


C° 


rl = -y-$U-M'U+v' x 


x ( 3 !+«?• 

With the given transformations the 
dovnwash for the inner space problem reads: 

X 


W K = — ^ " ■ * +Q <"- “*> [- \ («,. * 0 ) _ 


3 tj arc ctg cos # •e i,t . . (6,2) 3 . «*„<«., to f 

/ -~- nuaj-*r Q '(tux,) — — l _ 

i transformations the 


x 0 s • e 


— (oix »8 


ds 


+ 




w *-‘~ irr n ) x t 

4 V \ci i ^ ^ 

and that for the outer space: 

Co 

***.• = — f*’Y‘ j*( 3_ j-qp-^r — 3 ’? arcct K*?)e i,,r ‘<» 


+ 


cc 

l jlc 

* da >j 


x 0 s • e 


1 

-uix,M 


(** ** + 1 — *o) V «* 


(*?«*+ i_ x *)y s .-=rr 

.... (5,4) 


In this formula the integrals are 
already known from the calculation of 

Vi°. (See (4.5).) Moreover, this rep- 
resentation is beneficial for the com- 
pensation process later on. 

REPRESENTATION OF FLEXURAL OSCILLATIONS 


In the course of the symmetrical 
oscillations of the disk the oscilla— 
and with the substitutions given in the tions with parabolic deformation 
preceding section: 

GO 

= — yj*( 3 — 3 *»V** — 1 ;rc ctg V ** — 1 — 


W = D x2e iut and W = F y2 e iut , re- 


*o**+i — *0 




-wiXil, 




V**-x 


ds ■ 




_1_ f x, s ■ e~ “ i x ’ ‘ ds 

* 0 ** + ! — *o 1*~l' 


spectively, are dealt with. To pro- 
vide simpler expressions for the down- 
wash calculation, the following linear 
combinations are Instead used as basis: 
^*.1 = A •(*!+»*) I 
wVt = -Fi-(* , -y , )- « i,, ‘ I 
The individual deformations are 
obtained by additive superposition. 


(«. i) 
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Tho repre sentation of the first requisite dovnvash under ( 6 . 1 ) pro- 
ceeds from the pressure potential 


v . = C o£i!^(|^ + |_^^Larcctg,|. e ‘’-‘ 

The inner space problem yields the downvash 

X 

*lr- f*^F J [2 — 6(l'*+y«)] e ‘-‘' < ix' = 

— x. 

- Q \ Art l 3n , 16 jj , 15 n , . . „ . , _,\ 

~ V r ( TTeu HTw ITiaP 0,1 x * ^ 1 a> x-\- 2) J 

C *' [ ^4^+87^ y -87 ^ Ht0 ^ + 2 ltOX *+ 2 ) j/- 


5 n* — 3 ft [6 


2 5i7* + 3ij 


arc ctg 




and the outer space problem: 

I ??® f * 3 / s „*- 2 5 v * + 3 r/ 


"i— /5S*+S- 

0 

--Jill* 81 * 1 - 1 




— tai ( / x,* + »/* 

• 6 , d n = 

i *o + v 1 

(5 x\ (s 2 — 1) — 3) arc ctg a: 0 V s 2 — l] — 


Vs 2 — ! 


Then: 


That is: 


CD 

*0 j <s* — 1) (1 — *0 Vs*— 1 arc ctg x 0 V s 2 — l) yp=pds = — ( ^^ ° + ' ^i. • 

is: _ 

w t - = x (~3' + * w °- -) ~l W l°+i(~Y i ^ <“> *))• 


Insertion of the value according to (4.4) for fl?, a affords 


"1 - wi . + Vi. = % e irt -?-[-A + « w /^ + Si*_ a \ j 

V 4[ eoi 2iw ‘ 2io>\ co co 2 l\‘ 

I ^3 A(rt — IOX) ll , t 15*0 , 15 d* \ r n . , , , 

+T e l \- 1 + -i-+T d^][-T lH o (“*<>>] + 

CD 

, /_ 9 15ag 15 15 d 15 d 2 \ f x 0 se- aix ’ s ds \ 

T \ 4a>i 4 wi 2iw 3 2 ia> 2 d w 4 i co dm*! J a; 2 s 2 _j- i z 2 y s j j j * ■ ' (®>3) 

The first part of ( 6 . 3 ) yields the and of Bessel's differential equation 
required expression for the dovnvash of a* (2 . , 2 , i 

the bendina oscillation, combined vith d®' 0 OJ ' x ° ^ ‘ T ° ’ 0 m ' , a a,» h ° ■ ( a> ' *“) 


the bending oscillation, combined with 
a torsional and flapping oscillation, 
which shall now be compensated. 

By putting 


eg = w i *g 


changes ( 6 . 3 ) to 

w " = -e ir, + 


and forming + K l . e nn-,., x) / » . s wH w (eo . ^ _ 

IfO o j^O r \ZZ 

°i 5 f a: s-e~ iwx,s ds \ 

the utilization of the relation — v i .J* ,,, it ' < 6 - 6 ) 

CO J (*o s +1 — *o) Vs 2 — 1 / 

wi f • Vlfel = i n 2 i ~L + Formula ( 6 . 6 ) gives the requisite down- 
i ( 6 4 )Wash and an additive term to be compensated 

/i i d* \ Fx s-e~' oix ° s ds' later * second oscillation given under 

~ ~drt) \ 1 - a! ' T 7 =T 18 0,btain ® d -through the pressure po_ 

i -Mi tantia! 
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<rl = c\ (i + *!*) (is - - 


and for the outer space 


c 2 C 1 / 

— 15 >7 arc ctg >7) cos 2 £ • e 1 y ( (6,7) — -y J g- (1 + V*) (l5 

' <r 

The doimwash of the inner space then reads 


-2 


wi 


1 = ^jrJ — -jf -n [x^ — y 1 ) e ,,t ' dx' 

— X 3 


— 15 f? arcctg' 97 ] 


1 + iJ* (1+7*)* 
cos 2 ft j.c 


I 


dt], 


15 


ny 


-* e ivt 


W\ 1 = -=f [~ w t^-=( — C 0 7 x * — 2icax-(-2)-,- ltt . 

3 ‘ V [16 ico 3 ib ta> 

+ ?| e U- | -«'* ) fl| 3I .4- i -(+c«**S-2i«>*o-2) 

15 


vith 

COS ft : 


"cos2i- rrt ' +irt 


+ 


Vi-W 


1 + 7* 


- — 1 . 


■»r c 


lb i co 

or, after substitution in s: 


( 6 , 8 ) 


j s x 2 s 2 ‘e — tt} i x o s i f* - e — 

| (l5-15x 0 l's 2 — 1 arcctgXoVs 2 — l) - 0 y_—_ ds ~£ (1 — *0) J (l5-15x 0 Vs 2 — 1 arcctg a 0 V«*-l) yy ~ ds_ 


f/ 5 1 

2 \ (x 2 s 1 

J U*o + l-*o 

(xgs* + l-xg) 2 j 1 0 


— (otX,8 

- ds. 


With W?,_ and itB derivatives, the total downwaeh follows (cf. 
(6.2)) at: 

W\ = Wl, + Wl. 


/ 15 

16 (1 — *J) 

! 15 d 

15 d« \ | 

\4 i co* 

8 <0 i 

4 t <u* dm 

8 1 a) d <u* ' , 

J 

03 




1 (7 

6 1 

2 

\ 4* 

8 J Us 

;»* + i — *0 

(****+ 1 

— **)*/ 0 


x 0 s e 


r toix,8 ds 


— a > * X* « 


With (6.4) and ( 6 . 5 ) and the relation 

2 d C x„*-e- ,aix ’’ 

j y »*— 1 * » de °J (*Ss*+i— *S) v **— 1 


ft'H-l+l' *-«<*•» JV.-wi*.* 

1 ■ as = — 


+ l — *0 V** — 1 

it then affords 

CJ 15n j.| / 1 « . 1 . 2 ] 2 \ | Cj i( r |_o, x ) [/ 6 IS d IS d* \ / i im x )) 4- 


1 1» 

16 


1 15 d 

15 

_£_\ 

\ 8 co i 

4i«u» 

4i dco 

4 i co* <1 co 

8 i co 

dco*/ 


f 


+ 


Analogously 

*o»*— l + *o »' 


15 d» \ (* e -°» i3c -« ds 1_ f xgs» — 1 + xg e -joix,' 

ico d co*)J xgs* + l — *§V**— 1 4 J (xg** + l — xg) 2 Vs* — I' 


d s. 
• (6,9) 


■ xi*.* 


(*g»*+i — *g)* y^nrr 


■ <f # » 


-i[— / 


-ui (*«» + i Kl — *!) 


(*.» + < Vi-^)*V> 


For / = sin «; *, = cos a with allowance for 
00 


I ,, , n —o)i (*,« — » — 

=-d» + e + “f 1 -* : ( — 2 7 - ds ■ 

— l ^ J (*.s- iV 1 — xi) 2 Vs* — 1 


. ? — u> i (* • CO* a -j- i tin a) JP — to i {« • COS a +.» sin a) 

i -r- . d* = — i l — j== 

dm j (* • cos cx -j- * *‘ n “)* V ** — 1 J (» • cos a -f- i sin a) V s* — 1 

, cP — »i(i - cot a — iiino) /? 

_i_ \ - .. d, = -i \ — 

dm j)(*-cosa — ^sin^*)*|'* , — 1 J (» 


e — O) i (« • CO* a — i Sin a) 

• cos a — ism a) V ** — 1 


ds. 


ds. 
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It la 

? «S*-l+«» ,— «*« 
) (*S**+1— *o)* l'* 1 — ! 


, f . / J' C BJ(#-cot« + <«lno) \ . f f — a>U«-co«o— iaino) \T 

= i 1 C«>-i — : p=i,i5U. + ,, '"° C?'-i —5 -^=d*d<5 

2L \ 1 £ (#*coaa + i«m<*) V**— 1 / ' ^ g (# ■ cos a — isinac)V* a — 1 /J 


and, by similar considerations 


^ (*j^+i-xs* it-t aL i 1 1 V**— 1 J 


jjf 


- CO i (<5oa a • « — tsin a) 


Vs‘— 1 


5 (J 


ds d col (6,10) 


f- 


+ |C<*> — i C'j 11 • co — 

Constants c| 2 ^ and C^ 1 ) are again computed by changing to the 
s 2 

stationary value to >0 according to (4.6). 

These quantities posted in (6.10) followed by the integration in the 
triple integrals over a, which gives Hankel's function, result in 

ajP — 1 + sj a"”**.* 

(*|.*+ !—*«« 1/1^1 


- Oi (y) • ®of (a> sin a) + i — (y) • ©in {co sin «) — 




: (o> sin a) + <0 ■ Oo («/) • ©in {co sin a) + l?*®* (a> • COS a) -Eo[[(a> — co) sin a]da>. 


Putting 

Cf= icoAS 

and compensating the torsional oscillation still existing in (6.9) "by form— 
ing 

Wi = JVi + f~f-W!. 


the downwash formula (6.9) then becomes 

l(-T i “-Y i £-Ti < °-£t)(-I i H ?’{co*o)) 




+ 


. IS ,5 _d U5_dM f£o 

' 4 4 10 d co 4 lioi’/J ij 


ds 


'/ 


z»s a — 1 + sg e - wix ’ s 

(xgs' + l — *!)* y s a — 1 


ds 


»s 2 + 1 — is 2 — 1 4 

The term still appearing beside the requisite flerural deformation 
must be compensated later on. The downwash formulas are all changed to 
reduce the paper work to a minimum in the subsequent compensation process. 


( 6 , 11 ) 


The pure flexural oscillation in x or y direction could now 
equally be formed. But, since the downwash formula W* is very simple, 

while the formula for is much more complicated, it is advisable to 

postpone the linear superposition until after the compensation calculus. 


REPRESENTATION OF THE TORSIONAL OSCILLATION ABOUT THE X-AXIS 

(ROLL OSCILLATION) 

In connection with the oscillations symmetrical with respect to plane 
x - z, the unsymmetrical oscillations with respect to this plane with de- 
formation of the disk are analyzed. 
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The requisite downwash shall read: 

W — C' y-e ivt 

The calculation succeeds vith the pressure potential 

V« = £lyVl — >*-y — y^r — Sijarcctgn) sin^-e 4 "' 

The downwash formula of the inner space reads: 

__ * __ __ 

-J l*- JL e lrt + °*. e Hrt-»*).3 

• F 4 J F 4 co 1 ' F 4 011 ' 

— x. 

and that of the outer Bpace 


(7. 1) 
(7, 2) 






_5L<(W- 

~ F* 


: "““ x) ^( 3 ~ 3 *» " c ctg -^ir+V— 4) ‘ y ^r 


or, vith the previously defined abbreviations: 

n?>. - f » n . - i » J(4£fr=4, « 

Bearing in mind that 


ds 


with the constant 


T f do, 

^ *o s *+ 1 — x o V« s — 1 ^ *o s ’+ 1 — *0 V** — 1 


the total downwash follows at: 




F 


^ = w , S. 1 .’<+ w'gL 

2 3 y u( . C 2 i(v( — o>x) 1 3 zjt 2 ) 1 \ 

r T in +“F" e 


3y f ds , i .. f. f ds 1 „ , ,1 

20,^3 4 S 2 + 1_4 Vi^T + 2 j f) x 2 oS a + i-x§‘ y^nzi 2 Ql>(y, J 

It thus affords the requisite downwash (7.1) and a term to be com- 
pensated. 

BENDING OSCILLATION 

The last oscillation with a disk deformation of the second degree in 
x - and y shall be represented by the downwash 

W — E- x y-e irt 

It is solved with the pressure potential 

2 


?» = Cl - (1 + <J*) (is — 1 - _ | 5 v , - (i — 15 1? arc ctg qj sin 2# -e* 1,1 


.-P» -/*• + /« 

I - c, g 

The downwash formula of the inner space is: 

X 

w f . ’ = -y- j — -y- n x ' y • e ’ w **•' 

and with 

2 ^I + ^ Y 1 - 4 

l + *7* 


(7, 3) 


(8,1) 

( 8 , 2 ) 


sin 2 d = 
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that of the outer space reads: 


M-*o(**— 0 (*3**-f i-*S) a 


15 x» fs*— 1 are 


ctgxot'^ — lj 


XpS ■ t 


-wtx,* 


i *— i 


■ ds 


= T y 


^•+T«.fe^-7lfeT 


Also applicable is 


<2* 

l 

C *—•*•/ i 

1 \rf._ 

Y*-i 

4 i}/ 1 — xj 

y fe*— 1 \(x 0 s — t Vl— xj)* 

(XoS+jfl— X*)*/ 


i r 

CD 

{* to Sin a — CD i (tcosa — » Sin a) 

\ e e ds 

CO 

— cd sin a e — o> i (« COf m + i sin a) *1 


4 £ sin a [ 

J y ,1 j (* • cos a — i sin a)* 

J y ,« i ( s cos a + i sin a) 1 J 


which, after the usual differentiations and transformations, gives 

f *.*- e -“ 1 *** d, 

J (*8«*+ 1 — *S)* 

i 

fit cO M oo. ___ 

(T f ds • « - 5 8in « dS ) - (c<» _ i C<» a, - 
4 1 sin a \ 1 * — 1 / 4isina' J J J V s* — 1 > 

where and C^ s ) denote the values computed under (4.6). 

i l 

Hence 

00 

- idh (-<?.•©«"(“«"“)+ *• 2 <“ sina > + 

1 "> 

-f- Q t ■ m (Eof (a» sin a) — i — P a ©in (cu sin a) + i 1/ j, 21 (a) cos a) ©in [(<w — m) sin a] dm) . . (8,4) 


To obtain the requisite downwash put 

C,* = o>iJf 3 * 

form 


w% = w? + \w r B 


and obtain: 


[(*+£-£) (-1 


d s 


s 1 — 1 


* \ 

| f x 0 » • e * *• * ds . m i ( 

r x 0 ds i 

\4 4i dmj 

J xs»»+i-xgy,*_i 1 z 

\ (xg s * + i_*§)* y,,_ 1 J 


the requisite downwash with the term to be compensated 


Herewith all the downwash forms stipulated in the introduction are 
computed, and proceed to the potential functions of the second type to 
be used to effect the compensation of the superfluous terms of the down- 
wash function of the first type. 

Obviously other higher potential functions of the first type afford 
further disk deformations, whereby the downwash functions become continu- 
ously more complicated. 
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POTENTIAL FUNCTIONS OF THE SECOND TYPE 


On the other hand the third 
Cartesian coordinate x provides a 
very simple and comprehensive ex- 
pression as real part of a complex 
function 


These are obtained from those of the 
first type by a limiting process. For c d 'A* 1 _ . )n . t 1 \ . 

this purpose the linear combination Wn= ~YJz \ — 2 x'~ ( (»"+ iy) n (x'—i «)») g * " tdf> ' 

: ® J in k\ 

(9, 1) , 


„ _ »v!+i 


, &y;+i 


dc 


■S+(»+ «■ 


(9,5) 


is formed, where 
circular surface. 


a»j dc ' ' *' c '“’‘ with the relation 

c is the radius of the »’* + »* »* 

The also satisfy 1— /<* 


/** 


existing between x, y and 4. 

the potential equation. 

This differentiation could, of course, The downwash integral gives in 
be equally applied to the downwash of the 
potential functions of the first type for 
obtaining those of the second. But in 
this instance it is found to be simpler 
to compute the functions according to 
(9.1) and to define the related down- 
wash, because the final formulas follow 
a pattern. Again the pressure potential 

* iv( - const (9, 2) 

1b as in the 


<Pn = Vn-e 

is valid, whereby 4/ , 
stationary case, 

Ml -/»«)* 


the stationary only functions in y 
alone, that is, the order of magnitude 
of the upper limit plays no part; in 
this instance 4 can be disregarded 
with respect to unity. 

After this simplification of the 
downwash formula (9-5), the following 
substitution is made: 


u = j/l + - 


Vn — 


cos n # 
sin n# 


ft* cos’ « ’ 
where y = sin a 
On putting 


«0 


= / 1 


cos* a ’ 


ist. 


(». 3) _ , . 

gn (“) = ; 


Art 


(ft* + *?*) (1 + *?*) * sn (u cos a + £ sin a)" (u cos a — £ sin a) B ’ 

These potential functions approach downwash integral^ changes to 

infinity at the entire disk edge as 1 /n. c„ (—1)" d p g ( u ) 

"J W=F iu - 


Art 


The downwash then follows at 


iv _ C n . (-it" A. 
n ~ V 2 cos* « dz 


Wn = 


X d 
V dz 


z -»0 


I' 


Vnli", 2, t')dfl 


Formula (9.3) can then be written a 
gain in the coordinates u, y, z whereby 

cos nd mus t be developed in a series of 
sinnd 


(9, 4) after which the differentiation with re- 
spect to z and z >0 is readily car- 

ried out. It results in: 


_C„ 

V 


(-1) 


2 cos 1 


0" r f?n(“) — g n ( 1) j , 

[ J _ d “ -tnW+gn (OO) 
1 • • • (». 


Here g n ( 00 ) = o, 1 for n > 0. 

¥ith»(« — t') = a>[x — *') and 

g n = 

the downwash formula (9.6) becomes 

c« (-1) 


w* = 


V 2 cos’ 


+ s’)— /**»*] 2 — 
*• — 1 


1 


+ 


1 • • • (9, 7) 

Partial integration affords 


cos d. Formula (9.3) gives with cosnd 
Bind 

the symmetrical, and with sinnd the un— 
symmetrical, potential functions. For 
symmetrical functions the development 
then would give 

v« (/«, y, z) = c n {((i - 

- (") ft* V* [(1 - M*) (M* + **)-/*’ »*] 

+ (4V 4 y* [(1 ->*) (X** + *•)-/»* y’]^-- 2 + . .. ) j d “ = [-*"(“> • y u .“n] + 

, As this expression indicates, the 1 

prediction of the downwash would neces— _j_ 

sitate integration over the surface else 
the total expression for the operator 
_d 
dz 


dg„(u) 

du 


Vu* 


- d u. 


becomes Indeterminate . This 


=•2— >-0 

would lead to serious integration dif- 
ficulties in the nonstationary. 


Whence (9*7) becomes 

W - 1 . (-P"/ feft.(u) ^ du _ 

" 2 cos* oc \J du fu*—i 

1 

OO 

— £(*■'(») — gn (1)) • jrj=yj — gn (1) + gn M j. 
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that is, for any n>0: 
(-D' 


IV — J?®-. <»*) . 

* V 2 cos' 


With 


f “-^ <UC0,a - 1 (* iB aA , 

11" C dp u J ucosa — tsin a fu' — l cosa \2 0 / 

LL-l-ggg-. d „ ( 9 , 8 ) I _ - 

;*«.) iu J'u’-l t t Sintt f 

cosa J uc 


1. n — 0 . 


d u 


tcosa — t sin a J 1 u* — 1 


4 l5 - = 2ia.cos«e +<< " ucos “ 
a u 

y = sin a x' = w • cos a 
the downwash reads 

CD 


f „ / g ie - UC05 ° \ _ 

J \ (u cos a -f- i sin a)* / y' u* — 1 

1 Cj» + Cg» . since 
cos a 2 cos ex 2 


— iuu cos a 


_ (*< — »*) j* to) COS PC » g ** 

# V cos* a J y u* — 

and with Hankel's function 

GO 


-du (0,0) 


J 


A to U COS a 

-j ?-- dw — -j- -=- i ^ri l) («o • cos a) 
7 u* — 1 * 


The general calculation of the func- 
tions and is carried out by 

1 2 

the method employed in the preceding chap- 
ters and therefore not discussed in detail. 
It results in: 


(9.9) becomes 


c, «' 


,i (* t — ox) 


nr V P 

n * ~ V ■ cos*« 
Because 


<2 In . 




iJJj*' (co-cos 


«))• (9. 10) 


;«<”{*) = — 


G W = e o,sino( C (.)_|_ £ J , ( ^.- 7 i) e -»»in« 

0 

co 

G£° = e " sin 0 + 1> e + “ sin “dm) 


(9, 13) 


dx “» / 1 

the downwash formula finally reads where G v; denotes Henkel's function 

C e i(r, ~ •“*) JI C ) 1 

^ 0 = — jf — — icucosa-tf ' 11 (o> ■ cos «). ( 9 , 11 ) !1 . In this manner the individual 

2. »»1, 2 

With the abbreviating Btyle of writing G ^ B ^ cessi J el y built 

^ up with Hankel’s function and solved. 


ft. («) = 


A 10 u COS « 


A to U COS a 


( u cos a 4“ i ®> n <*/" ~ (** cos <* — * sin °0 n 

( e itOU cos a V 

(ucosa + isin a) B / 


The constants C 


(n) 


are defined 


the downwash reads according to ( 9 - 8 ): ^*+ 1 ) _ ± 

W, = 


according to formula (4.6). 
That means: 


cos a J \ (u cos a + i sin «)* 1 f u* — l 

• I ,io>u COS a . 

-ig-f.fi ( — — - — du (9, 

cos a J \ u cos a + item a 1.) u* — 1 


“ 81n “(-uMf P 2 *Ti<? 2 *+ 


12 ) 


and with 


Gi’ 


+ l ’“|Q 2 ft-i± » Y I>2k - 1 l 

2T“* (iT^zft — 2 ^ *^zft — 2 ) — 


, . 71 

2ft — 3 ± 1 y ^2k—3) ’ 


A to U COS a 


(« cos a + i sin a) n fu’ — 1 


w 

-f 

g,., = r ^ muco, ° 

2 J (» cos a — i sin a) n / «* — 1 

the integral can be expressed with 

f «-<* <UC0,a du 1 l* iam \ .. 

J ucosot+isina y 1 u* — 1 cosa \2 0 / 

isina f t 1 ”""- du 

cosa J u cos a + i sin a / u* — 1 


• = •••• + (— D* -i 2h+1 jj ■ • J^lH-{, 1, e+‘ usina dc5] 

0 — * * •* 

2fc-f 1 

G(ZA) = e ± WS ino ( _ 1} k J Qsifc _ i ± . » i , 2s _ j _ 

ltt> (t ^2 ft — 2 T ‘^2ft— 2 ) 

-2T a,, ( < ?2ft-3±‘Y^2ft- S )+ + 

+ (- 1)* • i 2k ff ... e ^ ra8fn » dco] (9, 14) 


2 A 



NACA TM No. 1096 


17 


The rest of the downvash formulas are built up with: 

; / .-<*co.. ' udu , cj— 1 ) +g£-» , gina C<’>-G< t -> 

\ (u cos a -)- i tin a)* / V u* — 1 C 0 «« 2 'cosa' 2' (9, 15) 

Whence the downvash (9.12) can be expressed as 

to ■ 

W * “ co?ir (*>• V ) -*<?.- Sin i ^ »(&)[[(«,— s) y] df5?-{- 

0. 

I / 6 * 

+ ‘ y ^ «°f (“ y) + Cl ■ ® f « (°> v) — * “ (y p » ®«" («o y) — »<?.• M («oy)j — ** J*J*^ »'«?* Sin [(a. — a>) y] d s} 1 — 

6 • 

_ M 

— ESSS [y * #o ' ~ * »» * {y Pf Sin (co y) — i <?. dot (a, y) + i JiL *J5T«> ©in [(o> — co) • y] d sj ] 


Now the relation (reference 3) 




(9,16) 


exists between the spherical functions P n which also holds for Q n . 
Herewith the above downvash formula becomes: 


nr » 


W l = 


2 dy 


(£of (coy) — i yy-- Sin (coy) + |i J*y go( [(<o — co) y]dm — 


Q> 

T UI ^ ) ® in [(“> — ") y]dW — iw^iH% ) +i*my fy iH™ Sin[(co — co) y]d5ll 
0 o ' 

Further transformation by writing 

to to to 

J*J> Sin [(co — co) y| d a> = coj* -Sin[(u> — 7S)y]dta — J*«> J^ 1 ' . ©in [(co — co) y] dco = 

0 # 0 o 

tO ft> 

= + a) J* ©in [(co — W) y] d CO + ■ I/O’— i P (jfjny + o>~ ) ®o[ [{«, — Jo) y] dm, 

0 y y^f \ oco / 

the application of the relation — Hq 1 ^ = _ cos a H ^ 1 ^ before (9.1l) 
affords the final formula ^ 1 

W'c = -y- (f 4^ ® of (eu y ) ~ ‘ ‘ © in ( " y) + -^£7 1 Y l ' ® < 5 ‘ 003 a > ®°f •[<" - “) yl «* s) • • • (9,17) 

The downvash reads: 


3. n = 2. 


^ 


C, e 


i(*t — cox) {“ 


V COS* o 


with 


W,. 


(' dg t [u) “ 

^ d« •y u ._ 1 du (9, 


18) 


dg t 

d i 


cos* c 


t,(u) = _ R I 2 cos <x . e ia>UCOSn \ i iw coscce* rauc °»° \ 
du \ (« cos a + i sin «)* / ■” \ (ucoscx -f- i sin «)*/’ 

fy [ — Qi®°f (co y) — iy p i ‘ @in (co y) -f- ‘ co |y i\i®of (co y) — iQ # Sin (co y)|+ 

to 

+ ** J*J* y * Hi, 11 ®»f [(co— co) y) d S -)- i sin a jy J>, Sin (co y) — iQ, ®of (co y) + 

+ ©' n (“y) + »y -Pi®nf (co y)j — yp co* (y P t Sin (coy) — iQ 0 £of (coy)j — • 

to 

~ £ *JXf T £ — “) v] dwj j + c *"- w |y P, ®o{ (cu y) — iQ„ ©in (co y) + 


0 * 


+ *J*y *'- H o 1, ®of [(co — co) y]dco + isin a|<?, Sin (coy) + i Pi®o[coy — 
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According to (9.1*0 and (9*15) the downwash (9.18) then reads: 

co 

— i co P a ©in (coy) — i Q„ Eof (co y)) — i* J*J y i ■ ©in [(m — m)y\dm j- j • 

The existing multiple integrals are reformed by repeated partial inte- 
gration as follows: 

h i?* h< 1)® 0 5 r(q> — cS) y] . 1 ft 5 'n-iiL,Mi) gtlI[(a> ~g) y ld5 

JJ ©in[(co — e5)y] d ' (n — 1)IJ ( ' 2 0 ©in [(to— co) y] 

o * 0 

n 

Conformably to (9.15) this can be consolidated in the downwash func- 


tion V in the following manner: 
n 

~(n — \yr J* -1 [(“> — “) 2 /] dc 


(n — 2)1 


J* (co — co) n -^o 1 ’ Eof [(to — co) jf] d co + 


(», 19 ) 


+ (" — ™) n ®in [(co — 5>) y]dw (n°—l) I S^° — 7S ) n ~ 1 ^ iIi o ) ©tn[(a>— ^ a>)y] d& 

o 0 

CO 



if the last two integrals on the left-hand side are combined and trans- 
formed through partial integration, by putting 

©in [(e> — to) y] dm = — — dffiof [(to — co) y] 

Further allowance for (9.16) and for the fact that 

H$? ) = — cos a H { } ) 
dm 0 1 

affords collectively: 

Wt - *> [U* • <So\ (»») + i y • ©in (co y) m '$2 l @i„ („ y) + i i ^ • Oof (to y)) + 

CO 

+ cos « S ^ y * ^ 1> ^ cos ^ ^ dZ5 \ (9,20) 

o 

xu 

The general construction of the n downwash formula is readily 
apparent from the foregoing. 

It is 

»F. = l)» + i [^/ - f i • Go( [(to-*) y] <*o + 

+ ^ 1 )] (f ' e °H«> V) — i ^ ■ ©i« (« V)) • n *"" 1 - <«* y) + < i f ^7 ' ©i" (« y)) •<*>"-* 

- (f ©in (• J/)) ( 17 ®°f (» ») + »' f - ® in ( “ y) ) • £un “‘+ •- 

.... (9,21 

These symmetrical downwash functions of the second type, which are 
even in y, are suitable for the compensation of the terms in the sym- 
metrical downwash functions of the first type (U.j), (5.*0, (6.6), and 

(6.11) which carry the factor e i( ut_an: ) ^ are for that reason unsuit- 
able for the representation of the airfoil motion. 


For the compensation of the corresponding terms of the unsymmetrical 
downwash functions of the second type (7.3), (8.5) the unsymmetrical down- 
wash functions of the second type which are odd in y are derived similarly. 


The potential unsymmetrical to the x - z plane reads 
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(9, 22) 


C * U * 1 + *?’) (1 + *?*)"'* Smn#e 
For computing the downwash the similar complex style of writing is 
resorted to 


W H = ^(- i) (-1)"^ 


i-*o 


1 f/ (1-^)" (!-/.«)» \ * iW ' . 

2 J \ (x' + i y) n (*' — i J7) n / *' ^ 


Application of the same transformations used in the prediction of 
the symmetrical downwash function, and putting 

l l 


<T n (u) = 


\(u • cos a -j- i sin a)" (ii-cosc* — i sin <x) n 
yields the corresponding downwash integral at 


) / _i Ct» u COS a \ 

= 2 / 7 - r-^—. — , 

\{u - cos a. * sin fx) n J 


y ( 2 COS* « J du 


(u) 


- du 


(9.23) 


Vu* — 1 

A 

The subsequent calculation of the integrals proceeds in the same 
manner as for the symmetrical downwash functions. The general term for 
the unsymmetrical downwash functions then reads: 

(At 

(a> (t)) n 1 n „ (ll 777\ a 777 _L- (9,24) 


f'n .»(>'< — M *; 

1 V e 


(— 1)« f - 

' [ cos a J 


(n 


^yp- a> 2 iH[ 1] ©in [(a) — a>) j/] d w + 


‘ (f Tjf ® in i (^=iyi «"" * (^7 ® n («■»)+* y ^ Kof («»))- 


(n-T)r \2 dy *> * dy (» — 2)1* 

17 6in * if Kof <" »))+ 7T=J)T <0 ”~‘(l7 ® in (<u ■ y) + i 2l7 £of (<oy) ) + • • • ] 


(» 


It is readily seen that this method allows a continuous said lucid 
presentation of the downwash functions of the second type, and further- 
more, that with cd — >- 0 the same stationary values of the downwash 
functions are obtained as in the cited report by Kinner. 


THE FLOW-OFF CONDITION (cf. reference l) 


The condition for smooth flow-off from the trailing edge implies 
the disappearance of the values of the potential functions for 

A =0 >/ = 0 — yC^X+g 

This condition is satisfied for every potential function of the 
first type, because they disappear on the entire disk edge. For those 
of the second type It had been stated previously that their values ap- 
proached infinity as l/u. 

For the linear superposition of the compensation process an infinite 
sum of potential functions of the second type is formed. Since this sum 
must approach zero at the trailing edge, the conditional equation is: 
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G py cos2k(p + ^ C 2 x+icos(2 X + l) <p = 0 
k=o X=o 

To evaluate this relation analytically, multiply (lO,l) "by 

cos( 2 X + 1 ) 9 , and integrate over cp from 0 to 5-. 

2 


( 10 , 1 ) 


Then 


4 . \ X+ ^ . . 

C 2 \+x = - (~1) (2X+1) 


(,-D k Csfc 


k=o 


(2 X + l) 2 — ( 2 k ) 2 


( 10 , 2 ) 


For unsymmetrical potential functions of the second type the condition 
of smooth flow-off reads: 


) Caksin 2 kcp + Y C 2 X+isin (2 X + l)<p = 0 (10,3) 

k=o \~o 

Multiplication of (10,3) hy sin 2 Kcp followed by integration gives 


Osk = — — (~l) k 2 k f 

n / 


(-1) C 2 X+i 

z - ( 2 X + l ) 2 - ( 2 k ) 2 
X=o 


CC’MPENSATION CALCULUS 


The general expression for the downwash function of the first type 
is according to ( 3 ,l): 


W . f (x,y) e Ivt + e l(vtr **> g(y) ( 11 , 1 ) 

where f(x, y) exactly represents the requisite downwash. The 

downwash functions of the second type are represented by 


W n = c i^n(y) 


e i(vt-u>x) 


( 11 , 2 ) 
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Equations (11,1) and (11,2) are linearly superposed, so that the follow- 
ing is complied with 


g(7) + Y C n h nfr) * 0 (11,3) 

whereby g(y) and h n (y) both must either be symmetrical or unsymmetrical 
downwash functions, respectively, the h^y) are generally read from (9,2) 
and (9,24). The odd C n (n * 2X + l) and even C n (n = 2 k), respectively, 
follow from the flow-off condition treated in the preceding section. For 
the stationary case (<o— > 0) the orthogonal characteristics of the 
spherical functions is employed in the analytical problem, resulting in 
infinitely many equations for infinitely many unknowns C n (reference l). 

In the nonstationary (o> ^ 0), no orthogonal functions to h n (y) 
are known. 

In view of the repeated occurrence of the terms in coe^a in (9, 21) 
and (9,24) the orthogonality of the trigonometric functions is employed. 

The problem could be solved for symmetrical g(y) and h n (y) (ll,3) 
by multiplying with cos 2ka and integrating from - 5. to + 

But has the general form 

dy 

+ „ lW mil 1 

dy l - y s 1 - y 

where ax(y), and bi(y) are polynomials in y. Since the denominator 
(1 - y a ) would impede the integration over a, multiply (11,3) by 

cos2a x cos 2k a and then integrate. It affords infinitely many inho- 
mogeneous linear equations for the still unknown C^. 

The unsymmetrical downwash functions of the first and second type 
are largely provided with the factor y * sin a but, for the rest, even 
again in y. So in order to enable the use of the integrations over a, 

multiply in this case (11,3) by - ? . ^ .. 2a cos 2k a and Integrate over a from 

sin a 

to + * , 

2 2 
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The result is again an infinite number of equations. for the still 
unknowns C 2 X+i • 

The constants C- and C , respectively, are, then, functions of 
to only. 

PRESSURE DISTRIBUTION, LIFT, AND MOMENTS 
From the relation 

b + i grad p » 0 (12, l) 

“ P 

imposed by Euler 1 s equation, vand b = grad cp follows the equation 

P 9 = Poo ~ P (12,2) 

where poo is the static pressure at infinity. With it the resultant 
pressure per unit of surface in a point at the airfoil for short pres- 
sure jump, is: 

II (x, y, t)Sp u - PqV = p (cp^ - q> u ) (12,3) 


where q = 0 in the potential functions. Since the values of the po- 
tential functions on the upper and lower surface are inversely equal, 
equation (12,3) can be written 

II (x, y, t) = p u -Vob = 2 P*Pob (12,4) 

Hence, with the values of the coefficients of the potential function 
defined by compensation calculus, the pressure distribution can be solved 
with (12,4) as it will be achieved numerically in part II. 

The lift follows at 

A = JJ* II Fd (12,5) 

with the surface integral extending over the disk of radius c. 
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With 

(12, 5) "becomes 


dF = c ududcp 

STT 1 

A s» q s / / XI pdjidcp 


( 12 , 6 ) 


p.= 0 

Bearing in mind the relation 


2TT 2 tt 

cos ncpd<p w 0 / sinncpdcp ss 0 n 


- 1> 2, 3 . . 


it is evident that 


Ai = | Ttpc 3 Ca e iVt A 0 ss 4rtp c 2 C p e ivt (l2,T) 

m 

while for all other ep n and cp n the lift disappears, so that the 

calculation of the total lift "becomes very simple. Thus the total lift 
will scarcely depend upon the convergence of the C n coefficients. 

This does not hold true for the pressure distribution. 

The pitching moment of the lift forces about the y-axis reads: 


2 TT 1 

M = j^j IT xdF = c 3 J J* II p. J 1 — p 2 coscpdpdcp (12,8) 

cp=o p=o 

With the relations given in the lift calculation, it affords in 
this instance 


m| as JL «pc 3 C 2 e iu ^ Mi = — «pc 3 Ci e* u ^ (12,9) 

15 3 

m 

while the moment for all other potential functions <p n and cp n vanishes. 

Naturally, this applies only to symmetrical potential functions. If 
they are unsymmetrical to the x, y-plane, the formulas (12,9) indicate 
the amount of the rolling moments, in which case the corresponding C of 
the unsymmetrical potential functions must "be inserted. 
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PART II 

By K. Krienes and Th. Schade 


This part deals with the aerodynamic coefficients for downwash dis- 
tributions which are arbitrary functions of the second degree, that is, 
in the range of the reduced frequencies a> = 0 to 2.0. The pressure dis- 
tribution is secured also in several cases. The numerical results are 
illustrated by graphical representation. 


INTRODUCTION 


As already explained in part I, the pressure or acceleration poten- 
tial is put down as sum of potential functions of the first and second 
type. The functions of the first type are obtained by a separation theo- 
rem from Laplace's differential equation = 0 after introduction of 
elliptic coordinates 


Vn (m, v, &) 


1 (n — m) I 

(n-f-m)! 


C?P?M Q?(ir,) 


cos m& irt 
sin m 0 


( 1 . 1 ) 


The 

of 

and 


potential functions of the second type are obtained by differentiations 
cp£[ +1 with respect to the base circle radius c for constants x, y. 


_ _ r A* (1 A* 1 ) 2 cosng irl 

” * sin n # 

(a**+*?*) U + n 1 ) 2 


( 1 . 2 ) 


The downwash velocities induced by thes& potential functions at the lift- 
ing surface — termed downwash for short — have been calculated in part I 
and are briefly restated here. 

The potential ?>, = £«,?>? (l. ®) 


00 

yields .r,— t-f *.«> + j 1, ] <>-) 


Vi= »“>?>* + 4- V 1 ? 


yields 


+ ■H* -3 A ) 1 4 *• + 1 =3 ' rfe ] 


9>1 = i o> vl — 5 v l — -S- i CO p? 


(1.5) 


(1.5) 

(1,7) 
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yields 


(The term to he compensated (cf. Compensation Calculation) is, up to the 
factor — ji, the same as for the doumrash function W g (1.4) 


4 — »u>4 + -§-4 


yields 


. . * a ig <** 1 r d» 

da, 4 <i<Wp*3»»+i-*s ‘y^rzr 

1. ( p ° 4 »*— 1+4 , i 

-4 ‘“J { 4*+i-& wl (1 - 10) 


of the potential functions antisymmetries! in y 


yields 




W t - - % ■ ± ■ y ’* + -Jr ( ” «■ “ - *> [| y - o, • f Hi" (® *.) 

3 , {** *.. «-”«*■ y. T ^“ <x - 1 

4 y**— 1 2 p 4 y«*— 1J 


4= ,a>4 + 4- 4 


_ A f g— « , < x .* <f, , /5 , i 16 d \ f yr* lM «1» 

4j4* , +i — 4 y«* — 1 \4 4 ido»/J 4* , + i — 4 y«* — 1 

yields m 


iA f •~“ tx ,§ d, 1 
' 2 p y^_ij 


Hereby 


f ^TT-T— j ’ 7 == = « <?• (V) 6in (“»>+T -P. (») *°f (• V) 

4 »*+i — 4 f* 1 — 1 2 

m ' 

— y — SJyJdS • (1. 1#) 

0 

«* 

— y io> .P, (y) ®of (roy) 4-°>C* (y) ®* n (°>y) +•][ *J*(® — ®) •®o*’ (® **) ®°f [(«» — ®) y]4® . . . . (1.16) 


2 6 
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1 (y) 


f «-«<«■» . f , f * *.«•«-”*** 1 n 

) *S*'+x-*5 ' y^-i + + y ^ 

00 

f d« — 1 (— <? t (y) ®in (g> y) 

P i*»** + X— *S)* jV-i *»y\ 

+ -y * P 1 (y) lof (o> y) +Q« (y) a> 4of (a> y) — y i P, (y) a> Silt (a> y) 

a> 

+ y i J* (a> — ®) [,*> (of a*) ® I" f (e» — B) y] d Usj 

0 

The symmetric dovnvash functions of the second type read: 

to 

W n = y - «* ( ” l “"* > (— I )« +1 [-^ J * 3 Y * M 1 ’ (" «of [(a) — S) y] i JET 

“--‘(f ■ («>!/)- i 4 f ® n { " y] ) 

- (Si “-‘(-tH *■»» v>) 

“ Si “ n -’(f ^ ~ f ^ ® in <“ y) ) + • • ] 

the antisymmetrical 

^ = 7 ^ (- i)" [-^J ( “ S’5r 55 f i *?’ <“ *»> ® in K® - s) y] d «r 

0 

+ j^=4)j “■■'(§ ^ ®» i» »> - * 4r &l <* “') 

- fSn' «-(4r *■ l "» l + t‘ '“"' l 

•] 


(X.I7) 


(1.18) 


U.l») 


(1.20) 


COMPENSATION CALCULATION 

The dovnvash functions of the first 
type can he generally vritten in the form 

w a = U (*, y) ‘ (y, o>) (w ~ “ *> . . (2,1) 

vhere f a (x, y) is a polynomial in x 
and y, while g a (y, cu) contains tran- 
scendental functions. This second term 
is "compensated" with the dovnvash func- 
tions of the second type, which are of 
the form 

W n = K(y,a>)e i * i -‘‘ “*> (2,2) 

The condition reads 

Ko,nW n = U(x,y)e i ' 1 . . . (2,3) 

n — 0 


that is, 

00 

— go (y,<o) = Z K <”" h n(y,<») • • • • (2,4) 
11=0 

To arrive at a numerical result for 
the coefficients K^, n, the infinite 

series is broken off with n = 4. Fur- 
thermore, equation (2.4) is multiplied 

by cos 2 a cos2kada (k = 0, 1, 2) in the 
case of the symmetric dovnvash functions 
and in the antisymmetric case by 

C ° S a cos2kada (k = O.l) and then inte- 
sin a 

grated from 0 to — . This affords a 

2 

first equation system for the K a , : 
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plex unknown is eliminated. Accordingly, e- 
quation (2.6) needs to be taken only for 


k = 0; 1. 


Again putting 


J*' 2 g„ (to, sin a.) C °* “ ■ cos 2 k<xda = g a k , 
0 

31 

J 2 h„ (to, sin *) ^ ~ cos 2k<xda = h n , k 
0 


(4,11) 


according to equation (3.2), equation 
(2.6) becomes 

1^1,* +*0,3 ^3.* = * — 0 ; a ( 4 , 13 ). 


with 


and 


^.4 — T“(— fr X -.» + T Z -*)i 


• (4,12) 


— — 8 — 16 — ' 
M \, k = *«, * + *«, » 

M i,k = ~ + 7 ^* 4 , k. 


. . (4. 14) 


The solution of equation (4.13) follows the same procedure as for the 
symmetrical case, except that all matrices are now of two rows. 


The next problem is to determine the functions g a ^, k , and so 

forth, ((4.2), (4 i 11 ) ) . To this end the downwash functions enumerated ear- 
lier are developed in series of < 13 , the coefficients being indicated as 
trigonometric series in a. 


The integration according to (4.2) and (4.11), yields complex series 
in ai for g_ , and h . , the values between o> = 0 and a> = 2 being 

(J jK 11 y iC 

computed in stages of 0.2. Then the real and the imaginary part of the 
equation coefficient M^jj- can be determined with the aid of (4.4) and 

(4.14), respectively. Since the solution of the coefficients of the o> 
series g 0 ^ and involves considerable paper work, it is well to 

know that the calculation of the integrals for the antisymmetrical func- 
tions is rendered substantially easier if the following relations are 
taken into consideration (cf. Introduction): 


cos 2 a n . d 

- cos 2 k (x — 
sin a aa> 


J (a> — cy)”" 1 n 

~(n — 1)1 


to — i H { i ] (a) cos a) (Silt [(eo — to) sin <x]deo 


cos 2 a 
sin a 


cos 2 k a J* | “ w 1 ^i 1> ( w cos <*) ©i n [( w — «>)■ sin a] d to 


and 


+ cos 2 <x cos 2 k <x W 

o 


m ■— i H { j* (to cos a) ©of [(to — a>) sin a] d a> (4, 15 


cos 2 x d 

|- — cos 2 ka. [f ( y ) ©in (a> sin a)] = cos* a cos 2 k txF (y) ©of (co sin at) (4, 16) 


In other words, the integrals of the antisymmetrical functions can be pre- 
dicted from the symmetrical functions by integrations with respect to o>. 

In the present report the relations (4.15) ancL (4.16) were used as check 
for the independently calculated symmetrical and antisymmetrical integrals. 
The results of the solution are tabulated and plotted. 
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DETERMINATION OF THE AIR FORCES ACTING Since only and cp contribute to the 

total lift (cf . part I), the lift during 
ON THE OSCILLATING AIRFOIL the flapping motion is: 


If the coefficients Ka, n derived in 

the foregoing section are entered in equa- 
tion (2.3), it results in 


A — n g V* c* • -jj- j^-i- o»* to — ito K' t o| • — ‘‘ ' 

The total moment, referred to the y-axis is 
given hy <p| and <Pi: 

M = n e V'<*- „ 8 [u> K‘. , — io> K\ J 

Zn 1 C 


W n +21 K«nW n ~f,( X ,y), i ''‘ . . . (5.1) 

n*= 0 

The single periodic motions of the circular 
airfoil given, say, in the form of 

= ( 5 . 2 ) 

can be computed then by linear combination 
of the f a (x, y). 

The Aownwash is: 


In the tables and curves the values of the 
dimensionless coefficients k A and m A 

are given conformably to Hissner *s definition: 

A =3i e V' c* • k. 

M = n e Po'-u. • — «<". 

c 

The calculations for the relaxation oscillation 
are as follows: 


*(.> - ^ = K [i m , n[x , y) + b ^ ] ) ’ 1 (5. 3) 

For the given flapping motion 
for example , , _ A e‘ ” 

C 

that is, lV A = V^- iwe iyt 


and by consideration of (l.h) and (5 • l) 
it is: 

Wj — V* • — i 0) ■ ft (x, y) e* ¥ 1 

A n c f 

— K* ■ ( a> I" W v y -Ka.n ^*1 

” c L n = 0 J 

Correspondingly, the respective' pressure 
potential is 

9a = — | y ■ [*+ 2! # **••"•] 

■= V«.A. iu> £ A.,nV»] 


The pressure Jump, that is, the lift den- 
sity at the surface is then given by 

n(x,y,t) = 2 e<Pct. 

which in this instance means 


n (x, y , «) 


4 k«.A. ®.L»y i-r» 

2 c si L 


— i 0) 


4 


2 

n= 0 


r-cosntf l iy< 

Vl-r» J ‘ 


whence, according to (5-3) 

W, = — — -V-(iwx + l)e iyt 

c 

= B • V' i <1 • A_ (h'* + „ IF.) 

+ B.V.±lw,+j£ K. tn W,\- 

' 71 = 0 ' 

With epjj. = iaxp| + | <P?, and ep 8 = ioxp? 
the related pressure potential reads 

+ 2" (xt. • +■3' * °* ^ »j v»j . 

from which pressure distribution, lift, 
and so forth, can be obtained again: 

A = n e v* «* • [— a:, 0 — Y “ **■ 0 

+ 0 + “ AI, 0) j ■ Be iyt 

+‘(A.i+-§-<» Ati)]- B-e i,t - 

The prediction of the flexural oscil- 
lations and of the airfoil motions anti- 
symnetrical in y follows the same pro- 
cedure . 
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*- f 

- ^ _ u . J i y l 

z k — e x y e 


W d = ~'Y • (»o>x’ + 2a:)e*' < ; 

W r = — ■Via>y*e i, ‘; 
c 

W c = -j-Vi(oye trt ; 
Ws^f-V-ttcox-y + y)'*' 1 . 


9c = - J • V* ■ ^ [t" 1 P 2 — i “> £ Kk. n 9n ] . - 

^ = f F, i[ - 2i “^+ T°**S 

n = l ' J 


The respective pressure potentials are: 


'* " T ' ^ ’A [(“ ! + 1 “•)' *! - 2 '■ • *4 - 
+±(-*..+£n.-£n.)4 



W = 7 - ’ V * ■ [— ®* 9% + y ®* ?? — 2 w* 

« 

+ i<o£ (AJ.. + 2JT?,.) y B . 

n = 0 J 


The potential functions cp n and cp n are 

given in part I. Their values at the cir- 
cular surface are obtained for T| = 0. The 
corresponding lift and moment coefficients 
are indicated with k D , kp, mp, and so 
forth. 


TABLES AND DIAGRAMS 


Tables 1 to 6 give the values of coefficients K q n divided in real 
and imaginary part, resulting from the compensation calculation in connec- 
tion with the flov-off condition. The imaginary part K* is zero for 

„ j- o.n 

co = 0. ’ 


The framed— in values are equal to 



The real and Imaginary part of the lift and moment coefficients are 
also given for the different oscillation modes; the graphical representa- 
tions show these in complex representation in the Gaussian plane of num- 
bers . 


Figures 11 to 20 show in perspective representation the real and im- 
aginary part of the pressure at co = 0, 1.0, 1.8 referred to dynamic pres- 
sure - V 2 and for the amplitudes — anv - = 1, and so forth, respec- 


2 c o 

tively. The values are plotted for the wing sections 




= 0.383; = sin — = O.7O7; sin = 0.924 of the right and left half of the 

4 8 

circular airfoil. Incidentally, it should be noted that the pressure dis- 
tribution at the trailing edge of the airfoil was extrapolated with respect 
to zero, which does not follow mathematically because of the finite number 
of the series terms (cf. Kinner, reference l). 
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Figures 21 to 24 give the comparison of 'the presented calculations 
with the values frcm the strip theory- and -from Possio’s theory (reference 
6 ). By the strip theory each wing section is computed according to the 
two-dimensional theory without regard to the mutual interference of the 
sections and in particular to the finite span. According to it the sta— 
dc a 

tionary value iB = 2;t « $.28; according to Possio, who obtains the 

da 


dc ft 2 itl 

formulas of the lifting line for e» = 0 , it is — — = ~ 2.44. The 

da 1+2 


premises for his calculation were: great aspect ratio X and small to 
values. According to circular airfoil theory (reference l) the stationary 
dca „ 

value is — — = 1 . 82 . 
da 


Translation "by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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Table I. Flapping oscillation: ^ = 
Coefficients k s ,„ = r s , n + i k‘,. 


to 

0 

03 

0,4 

0,6 

03 

1.0 

13 

1<4 

1,6 

1,8 

2,0 


0^631 

03456 

03261 

03905 

03413 

0 t 1786 

0,1056 

0,0276 

— 0,0484 

— 0,1139 

— 0,1576 


—0,5489 

— 03367 

— 0,4988 

— 0,4367 

— 03487 

— 03373 

— 0,1086 

0,0281 

0,1613 

03789 

0,3686 

*;* 

03246 

03155 

0,1891 

0,1450 

0,0837 

0,0086 

— 0,0753 

— 0,1609 

— 03412 

— 031 16 

— 0,3761 


0,0042 

0,0090 

0,0234 

0,0466 

0,0766 

0,1098 

0,1416 

0,1668 

0,1806 

0,1797 

0,1640 

*1,4 

— 0,0475 

— -0,0493 

— 0,0547 

— 0,0624 

— 0,0699 

— 0,0743 

—0,0720 

— 0,0690 

— 0,0308 

0,0183 

0,1020 


t. S. 288 
— 03484 

— 0,0622 

— 0,1083 

— 0,1646 

— 03164 

— 03563 

— 03787 

— 03813 

— 03604 

— 03151 

— 0,1462 

*i» 

0,4465 

0,0931 

0,1913 

03879 

0,3737 

0,4387 

0,4739 

0,4725 

0,4314 

0,3606 

03330 


— 03062 

— 0,0622 

— 0,1243 

— 0,1813 

— 03271 

— 03566 

— 03617 

— 03432 

— 03006 

— 0,1370 

— 0,0682 

*1. 

0,1236 

0,0241 

0,0448 

0,0587 

0,0626 

0,0639 

0;0311 

— 0,0052 

— 0,0524 

— 0,1069 

— 0,1599 


0,0078 

0.0023 

0,0078 

0,0183 

0,0356 

0,0604 

0,0927 

0,1312 

0,1740 

03186 

03607 


Lift 

and moment coefficient <*4 + m A 

= m'j + im '/.) 




0 

0,0073 

0,0255 

0,0543 

0,1044 

0,1986 

0,3705 

0,6608 

1,1120 

1,7641 

2,6607 


a 

—0,1760 

— 03311 

— 0,4439 

— 0,4916 

— 0,4647 

— 0.3226 

— 0,0984 

0,1970 

0,6222 

0,8025 

m’ A 

0 

0,0158 

0,0660 

0.1466 

03638 

0,3724 

0,4827 

0,5615 

03859 

0,5357 

0,3956 

**A 

0 

03909 

0,1694 

03224 

03368 

03014 

0,1106 

— 0,0334 

— 03191 

— 0,4261 

— 0,6255 


Table II. Relaxation oscillation: = — 

C 

Coefficients jr J5 .-r«. + iii ; , 


to 

0 

03 

0,4 

0,6 

03 

13 

13 

1,4 

13 

13 

2,0 

^r.O 

—03221 

— 03212 

— 03133 

— 0,1941 

— 0,1613 

— 0,1153 

— 0,0683 

0,0049 

0,0681 

0,1247 

0,1661 

*r.l 

03872 

03844 

03677 

03299 

03681 

0,1831 

0,0797 

— 0,0336 

— 0,1462 

— 03470 

— 03258 

AS-., 

—03446 

—03397 

— 03213 

— 0,1861 

— 0,1338 

— 0,0668 

0,0100 

0,0895 

0,1646 

03282 

03824 

*K.t 

0,0888 

0,0829 

0,0648 

0,0351 

— 03041 

— 0,0490 

— 0,0946 

— 0,1345 

— 0,1636 

— 0,1767 

— 0,1729 

*X,4 

0,0069 

0,0117 

03252 

0,0453 

0,0694 

0,0937 

0,1138 

0,1251 

0,1223 

0,1013 

0,0499 

n, 0 

8. S. 28# 
0,1259 

0,0302 

03686 

0,1113 

0,1626 

0,1865 

03074 

03107 

0,1934 

0,1643 

0,0927 

Air.i 

—03830 

— 0,0603 

— 0,1330 

—03109 

— 03844 

— 03433 

— 03780 

— 03814 

— 03499 

— 03839 

— 0,1867 


03697 

0,0549 

0,1132 

0,1698 

03181 

03516 

03648 

03660 

03224 

0,1724 

0,1124 

Ajr,t 

— 0,1934 

— 0,0382 

— 0,0730 

— 0,1006 

— 0,1171 

— 0,1191 

— 0,1045 

-0,0734 

— 0,0283 

0,0266 

0,0819 

^,4 

0,0887 

0,0165 

0,0286 

0,0332 

0,0280 

0,0111 

— 0,0179 

— 03586 

- 0,1092 

— 0,1683 

— 03354 




Lift and 

moment coefficients 






03992 

03699 

0,7813 

0,6264 

0,4073 

0,1381 

— 0,1537 

— 0,4304 

— 0,6484 

— 0,7617 

— 0,7160 


0 

0,1314 

03684 

0,4020 

03904 

03517 

13085 

1,6720 

23380 

23889 

33835 


— 0,4659 

— 0,4524 

— 03114 

— 03398 

— 03397 

— 0,1203 

0,0015 

0,1043 

0,1640 

0,1592 

0,0703 

mi' 

0 

0,1225 

03466 

(V3664 

0,4386 

0,4760 

0,4563 

03746 

03339 

0,0460 

— 0,1710 
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Table III. Flexural oscillation 2: ^5=^ (**—»* )«**'. 

Coefficients Jti , n = x *!'!+» xl™- 


r co 

0 

03 

0,4 

0,6 

03 


13 

1,4 

1,6 

13 

2,0 

AT ? 

0,1123 

0,1077 

03937 

0,0701 

0,0372 

— 0,0034 

— 0,0487 

— 0,0949 

— 0,1378 

— 0,1767 

— 03090 


— 0,2620 

— 03531 

— 03286 

— 0,1823 

— 0,1212 

— 0,0471 

0,0338 

0,1133 

0,1825 

03349 

03655 

■ ns ? 

0,3034 

03944 

0383 

0,2257 

0,1685 

0,1015 

0,0311 

— 0,0341 

— 0,0860 

— 0,1121 

— 0.1036 

-nv 

— 03465 

— 0,2369 

— 0375 

— 0,1602 

— 0,0982 

— 0,0256 

03616 

0,1269 

0,1936 

03466 

03764 


0,1144 

0,1067 

0,0776 

0.0322 

— 0,0271 

— 0.0978 

— 0,1759 

— 03586 

— 0,3429 

— 0,4298 

— 0,6243 

* V 2 

S . S. 288 

— 0,1580 

— 0,0323 

— 0,0646 

— 0,0946 

— 0,1185 

— 0,1334 

— 0,1366 

— 0 , 1283 , 

— 0,1030 

— 03678 

— 0,0217 

Jsrj ' i ' 

0,3199 

0,0648 

0,1286 

0,1863 

0,2313 

0,2676 

03602 

03369 

0,1884 

0,1202 

0,0402 

JfS‘2 

— 03360 

— 0,0870 

— 0,1305 

— 0,1851 

— 0363 

— 03466 

— 03425 

— 03147 

— 0,1634 

— 0,0979 

— 0,0323 

*V,° 

03433 

0,0678 

0,1304 

0,1826 

03197 

03374 

03334 

03069 

0,1592 

03956 

0,0229 

Jt”2 (•) 
A B, 4 

— 0,2817 

- 0 , 0656 - 

— 0.1066 

— 0,1489 

— 0,1794 

— 0,1950 

— 0,1944 

— 0,1768 

— 0,1430 

— 0,0909 

— 0.0137 


Bending oscillation 1: wl=D 1 (x t +y*) « irt - 
For the coefficients/TA*the relation (c.f (i. 8)): *i. „ = — 
is applicable. 


Table IV. Flexural oscillation: z D = ^ ; x 2 e irl . 

Lift and moment coefficients 


“ 

0 

03 

0,4 

0,6 

0,8 

1,0 

13 

1.4 

1.8 

13 

2,0 


— 0,9436 

— 03631 

— 1,0005 

— 1,0999 

— 13619 

— 1,4439 

— 1,6603 

— 1,8369 

— 1,9598 

— 1,9824 

— 1,8568 


0 

0,0291 

0,0970 

03012 

03340 

0,4840 

0,6363 

0,7727 

03720 

0,9129 

0,8661 

m D 

— 0,4382 

— 0,4268 

— 03828 

— 03992 

— 0,1770 

— 0,0248 

0,1414 

0,3002 

0,4289 

0,5071 

0,5188 

Tflj ) 

0 

03202 

0,0464 

0,0660 

0,0606 

0,0136 

— 0,0878 

— 03491 

— 0,4675 

— 0,7312 

— 1,0239 


Bending oscillation: *,= — • 

C 

Lift and moment coefficients 


ft ) 

0 

03 

0,4 

0,6 

03 

1.0 

13 

1.4 

1,6 

13 

2,0 

K 

0 

0,0023 

0,0079 

0,0158 

0,0207 

0,0436 

0,0718 

0,1185 

0,1928 

03042 

0,4603 

K’ 

0 

— 0,0441 

— 0,0849 

— 0,1194 

— 0,1437 

— 0,1539 

— 0,1472 

— 0,1230 

— 0,0841 

— 0,0407 

— 0,0163 

m' F 

0 

0,0023 

0,0100 

0,0236 

0,0427 

0,0668 

0,0902 

0,1121 

0,1271 

0,1296 

0,1137 

m F 

0 

0,0189 

0,0359 

0,0494 

0,0670 

03606 

0,0461 

0,0248 

— 0,0070 

— 0,0463 

— 0,0883 
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Table V. Roll oscillation: z c = ^-yj". 
Coefficients k „, . = X;,» + i T k , .. 


CO 

0 

0,2 

0,4 

0,6 

0,8 

1,0 

1.2 

1,4 

1,6 

13 

2,0 

K ‘ s . i 

0,2264 

03177 

0,1065 

0,1625 

0,1197 

0,0734 

0,0300 

— 0,0045 

— 0,0256 

— 0,0324 

— 0,0254 

*«. 2 

— 0,2448 

—03351 

— 0,2074 

— 0,1638 

— 0,1090 

— 0,0497 

0,0067 

0,0627 

0,0828 

0,0940 

0,0889 

*«.3 

0,1050 

0,0987 

0,0798 

0,0508 

0,0146 

— 0,0249 

— 0,0631 

— 0,0960 

— 0,1201 

— 0,1325 

— 0,1324 

**,4 

0,0001 

0,0021 

0,0086 

0,0182 

0,0300 

0,0430 

0,0561 

0,0683 

0,0787 

0,0854 

0,0877 

^*.1 

S . S . 286 
— 0,1876 

— 0,0368 

— 0,0702 

— 0,0962 

— 0,1109 

— 04118 

— 0,0981 

— 0,0721 

— 0,0390 

— 0,0029 

0,0294 


03593 

0,0509 

0,0968 

0,1322 

0,1525 

0,1546 

0,1374 

0,1039 

0,0594 

0,0094 

— 0,0383 


— 0,1965 

— 0,0387 

— 0,0730 

— 0,0993 

— 0,1146 

— 0,1170 

— 0,1063 

— 0,0837 

— 0,0516 

— 0,0136 

0,0263 


0,0793 

0,0157 

0,0292 

0,0396 

0,0457 

0,0471 

0,0440 

0,0364 

0,0243 

0 , 0089 . 

— 0,0091 


Rolling moment coefficient m c — me 4 - i m'c- 


m c 

0 

— 0,0004 

— 0,0022 

— 0,0081 

— 0,0222 

— 0,0499 

— 0,0963 

— 0,1647 

— 03544 

— 0,3638 

— 0,4859 

m 

0 

0.0246 

0.0445 

0,0562 

0.0542 

0,0416 

0,0204 

— 0.0036 

— 0,0231 

— 0,0330 

— 0,0287 


Table VI. Torsional oscillation: z E = ~-- x - ye xrl - 


Coefficients a%„= kV? + i L‘. 


CO 

0 

03 

0,4 

0,6 

0,8 

1,0 

13 

1,4 

1,6 

1,8 

2,0 


— 03431 

— 0,2328 

— 03069 

— 0,1623 

— 0,1068 

— 0,0467 

0,0093 

0,0526 

0,0758 

0,0766 

0,0553 

/ r 3 W) 

^ B, 2 

0,3301 

0,3148 

03721 

03045 

0,1190 

0,0268 

— 0,0629 

— 0,1356 

— 0,1821 

— 0,1973 

— 0,1805 

V 3 (') 
■ n -«. 3 

— 0,2430 

— 03301 

— 0,1911 

— 0,1310 

— 0,0556 

0,0270 

0,1080 

0,1785 

03312 

0,2598 

03622 

QV 

0,0943 

0,0884 

0,0691 

0,0402 

0,0042 

— 0,0355 

— 0,0754 

— 0,1120 

— 0,1425 

— 0,1631 

— 0,1720 


s . S . 286 
03344 

0,0458 

0,0873 

0,1187 

0,1346 

0,1314 

0,1063 

0,0650 

0,0120 

— 0,0456 

— 0,0966 


— 0,4101 

— 0,0804 

— 0,1528 

— 03084 

— 03396 

— 03410 

— 03101 

— 0,1526 

— 0,0747 

0,0132 

0,0979 

xV ? 

0,4145 

0,0817 

0,1545 

03114 

03460 

03542 

03354 

0,1913 

0,1267 

0,0500 

— 0,0313 

xl ! -2 

— 03220 

— 0,0439 

— 0,0828 

— 0,1135 

— 0,1333 

— 0,1403 

— 0,1352 

— 0,1171 

— 0,0881 

— 0,0519 

— 0,0100 


Rolling 

moment 

coefficient mE = m e + i m ' E ' . 





m E 

0,1276 

0,1211 

0,1033 

0,0769 

0,0434 

0,0118 

— 0,0119 

— 0,0232 

— 0,0188 

0,0002 

0,0294 

m' E 

0 

0,0139 

0,0322 

0,0593 

0,0990 

0,1525 

03186 

03927 

0,3676 

0,4370 

0,4944 


34 
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Figs. 1,2, 3, 4, 5, 6 





f 












Figure 1.- fi and 5 indicate -tts — — 1 1 1 1 1 

P ositiv ® Figure 2.- Flapping oscillation A 
direction of the moments. ° lift coefficients in 

complex representation. 

I *i /<*•• I I 



Figure 3.- Flapping oscill- 
ation A; moment 
coefficients in complex rep- 
resentation. 


Figure 4.- Relaxation os- 
cillation B-x, 
lift coefficients in com- 
plex representation. 

jb'oij'W 



« 

" -t* -te •%* -it - to o 

Figure 6.- Bending oscillation 
D-x 2; lift coeffi- 
, cients in complex representation 

Figure 5.- Relaxation oscill- 
at ion in B-x, mo- 
ment coefficients in complex 
J representation. 
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Figs. 7,8,9,10,11,12,13 
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plex representation. 



Figure 8.- Bending oscill- 
ation F*y 2 ; lift 
coefficients in complex 
representation. 

02 Figure 10.- 
Rolling 
oscillation 

moment in 
l<t complex rep- 
resent ation. 


I '"J— i 

-63 -a 4 -qi <2 
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C-y; rolling 

VMAMM 4 
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Figure 9.- Bending oscill- 
ation F*y 2 ; mo- 
ment coefficients in com- 
plex representation. 




Figure 12.— Flapping oscillation 

A; limit case of im- real pree 
aginary pressure distribution. at w = l. 
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Figs. 14,15,16,17,18,19 



Figure 14.- Flapping 
oscilla- 
tion A; imaginary- 
pressure distribu- 
tion at w = 1. 





Figure 16.- flapping 
oscilla- 
tion A; imaginary 

j - 2 pressure distribution 

Figure 15.- Flapping w = 1-8. 

oscilla- 
tion A; real pressure 
distribution at tu = 1.8. 


Figure 18.- Relaxa- 
tion 

oscillation B*xj real 
pressure distribution 

at w = 1.0. 


t-ion 

oscillation B-x; real 
pressure distribution 
at co = 0. 



Figure 19.- Relaxa- 
tion 

oscillation B-xj 
imaginary pressure 
distribution at 

u> = l.o. -*• 



I 
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Figs. 20,21,22,23,24,25 



Figure 21.- 
Relaxation 
oscillation 
B-x } imagin- 
ary pressure 
distribution 
at to = 1.8. 


Figure 20.- Relaxation 

oscillation B-x; 
real pressure distribution 
at tu * 1.8. 




o ~—x 


-Figure 22.- Flapping oscillation 
A; lift coefficient 
k A i (real part) plotted against 

present calculation, 

Possio's theory, 

strip theory (two-dimen- 
sional flow). 

<«r 


Figure 23.- Flapping oscillations- -» 
A; lift coefficient 
k A n (imaginary part), (cf fig. 21.) 
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^Figure 24.- Relaxation oscilla- 
tion B-x: lift coeffi- 
cient kfi* (real part). (of fig. 21. 


Figure 25*- Relaxation oscilla- " 
tion B-X} lift coeffi- 
cient Icb" (imaginary part). 





